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Abstract. We examine the behavior of the coefficients of powers of polynomials 
^f) over a finite field of prime order. Extending the work of AUouche-Berthe, 1997, we study 

f^*) o.{n), the number of occurring strings of length n among coefficients of any power of 

a polynomial / reduced modulo a prime p. The sequence of line complexity a{n) is p- 
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regular in the sense of AUouche-Shalit. For / = 1 + a; and general p, we derive a recursion 
relation for a{n) then find a new formula for the generating function for a{n). We use 
the generating function to compute the asymptotics of a{n)/n'^ as n — >■ oo, which is an 
explicitly computable piecewise quadratic in x with n ~ [p"^ /x\ and a; is a real number 
Q^ between 1/p and 1. Analyzing other cases, we form a conjecture about the generating 

1 -p, function for general a{n). We examine the matrix B associated with / and p used to 

compute the count of a coefficient, which applies to the theory of linear cellular automata 
'"^ and fractals. For p ~ 2 and polynomials of small degree we compute the largest positive 

'^ eigenvalue. A, of S, related to the fractal dimension d of the corresponding fractal by 

^^ d — log2(A). We find proofs and make a number of conjectures for some bounds on A and 

^^ upper bounds on its degree. 

^' 
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^^ 1 Introduction 



X 



It was shown by S. Wolfram and others in 1980s that l-dimensional linear cellular automata lead 
at large scale to interesting examples of fractals. A basic example is the automaton associated to a 
polynomial / over Z/p, whose transition matrix Tj is the matrix of multiplication by /(x) on the 
space of Laurent polynomials in x. li f = 1 + x, then starting with the initial state go^x) = 1, one 
recovers Pascal's triangle mod p. For p = 2, at large scale, it produces the Sierpinski triangle shown 
in Figure M Similarly, the case of/ = 1 + x + x^, p = 2, and initial state go{x) = 1 produces the 
fractal shown in Figure [2] 



The double sequences produced by such automata, i.e., the sequences encoding the coefficients 
of the powers of /, have a very interesting structure. Namely, if p is a prime, they are p-automatic 
sequences in the sense of [3j. In the case / = 1+x, this follows from Lucas' theorem that (^) = Yl (^0 

i 

mod p, where rii, ki are the p-ary digits of n, k. 

In [6l[7], S. Wilson studied this example in the case where / is any polynomial, and computed the 
fractal dimension of the corresponding fractal. The answer is /3 = log (A), where p < X < p^ \s the 
largest (Perron-Frobenius) eigenvalue of a certain integer matrix B associated to / (in particular, 
an algebraic integer). In terms of coefficients of powers of /, this number characterizes the rate of 
growth of the total number of nonzero coefficients in /* for < i < p": this number behaves like 
n^ . The number of nonzero coefficients of each kind can actually be computed exactly at every step 
of the recursion, by using a matrix method similar to Wilson's; this is explained in the paper [3]. 

In this paper, we compute the eigenvalues A and their degrees for p = 2 for Laurent polynomials 
/ of small degrees, observe some patterns, and make a number of conjectures (in particular, that 



A can be arbitrarily close to 4) in Section 3.3 We also prove an upper bound for A depending on 
the degree of /. 

The size of the matrix B (which is an upper bound for the degree of A) is the number of accessible 
blocks (i.e., strings that occur in the sequence of coefficients of /* for some i) of length deg{f) (for 
p = 2). This raises the question of finding the number a(n) of accessible blocks of any length n. The 
number a{n) characterizes the so-called line complexity of the corresponding linear automaton, and 
is studied in the paper [T]. It is shown in [I], [5], and references therein that Ci'n? < a{n) < C2n^, 
and that for p = 2 and / = 1+x, one has a{n) = in? — n + 2. More generally, however, the sequence 
a{n) does not have such a simple form, even for / = 1 + x and p > 2. The paper ^J derives a 



recursion for this sequence, and we derive another one in Section 2.2. 1[ which is equivalent. These 
recursions show that the sequence a{n) is p-regular in the sense of |^ (the notion of p-regularity is 
a generalization of the notion of p-automaticity, to the case of integer, rather than mod p, values). 



We then proceed to find a new formula for the generating function for a{n) in Section 2.3, and use 
it to compute the asymptotics of a{n)lv? as n — ^ oo in Section 2.4 It turns out that if n = [p^/xj , 



where x is a real number between Xjp and 1, then f{n)/ii? tends to an explicit function of x, which 



is piecewise quadratic (a gluing together of 3 quadratic functions, which we expHcitly compute). In 



Section 2.4 we also compute the maximum and minimum value of this function, which gives the 
best asymptotic values for Ci and C2. This gives us new precise results about the complexity of 
the Pascal triangle mod p. We also perform a similar analysis for f = 1 + x + x'^ and p = 2, and 
make a conjecture about the general case. 




Figure 1: Fractal corresponding to 1 + x modulo 2 (Sierpinski's Triangle) 










Figure 2: Fractal corresponding to 1 + x + x^ modulo 2 



2 Accessible Blocks 

2.1 Definitions 

A block is a string of mod p digits. An ?7i-block is a block with m digits. For example, the four 
2-blocks modulo 2 are 00,01,11, and 11. 

For a polynomial /(x) with integer coefficients reduced modulo p, an accessible m-block is an 
m-block that appears anywhere among the coefficients, ordered by powers of x, of powers of f{x) 
modulo p. The number of accessible 0- blocks we define to be 1. Furthermore, we define row k for 
some f{x) and p to be the coefficients of f{x) reduced modulo p and define ajuYp^m) to be the 
number of accessible m-blocks for the polynomial /(x) and prime p. 

Example 2.1. For f{x) = 1 + x and p = 2, the 4-blocks 1101 and 1011 are never a substring of 
any power of 1 + x reduced modulo 2. Every other 4-block appears in some power of 1 + x reduced 
modulo 2, so ai+2,^2(4) = 14. 

2.2 Recursion Relations for a{n) 

We start with the well known fact in Lemma l2.2[ 
Lemma 2.2. f{x)^''P = f{xP)^ (mod p). 



Applying Lemma 2.2 to the accessible blocks, we have Corollary 2.3 



Corollary 2.3. For any integer k, prime p, and polynomial f{x), every row k ■ p for f{x) mod 

p is of the form 6iO . . . O62O 06„_iO . . . 05^ where the entries hi are the coefficients of f{x)^, 

and where each string of zeros between two entries hi and 6j+i is of length p — 1. Therefore, every 
accessible block from a row divisible by p is a subsection of 61O . . . O62O 06„_iO . . . 0&„. 

2.2.1 Accessible m-Blocks for f{x) = 1 + x and General Prime p 

The number of accessible TTi-blocks for /(x) = 1 + x and any prime p, ai^x,p, is defined by the 
recurrence relation in Theorem 12. 4[ 



Theorem 2.4. For f{x) = 1 + x and any prime p >2>, for < /c < p — 1, the recursion relation 
with starting points ai+a;,p(0) = 1, ai+x,p(l) = p, and ai+a;,p(2) = p^ is 

I ,\ {p — kMp — k ^ \) / \ /, , ,2 P^~Pn / 1N 

ai+x,p(p • n + A:) = ai+a;,p(n) + (fcp + /c - fc H — ) • ax^x,p\n + 1) 

+ -^!— ■ ai+a;,p(n + 2) - (2p - \){2p - 2). 



Proof. From Corollary 2.3 every accessible block in a row r with r = (mod p) is formed by 
adding p — 1 zeros between every digit of an accessible block, then adding some number of zeros 
(possibly none) less than p to either side. Furthermore, because f{x) = 1 + x, the coefficient of x* in 
a row is the sum modulo p of the coefficients of x* and x*~^ in the previous row. Because accessible 
blocks are subsections of a row, any accessible m-block comes from an accessible (m + l)-block. 
Table [I] provides the general forms of the {p ■ n + k)-hlocks for each row modulo p. To count the 
multiple additions of b in the forms, we define gi = (^~ ). 

The number of accessible blocks that lead into each form in Table [T] are the triangular numbers 
counting downwards for ai+x,p{n), the triangular numbers counting upward for ai^x,2{n + 2), and 
because the total number of forms is p^, we find ai_|_a;^p(n + 1) through subtraction. Namely, the 
factor of ai^x,p{'n) starts at p for row congruent to modulo p and k=0, and decreases as k and 
row increase, and the coefficient of ai+2;,p(^ + 2) starts at for row congruent to and 1 modulo 
p and increases with k and row. An additional (2p — l)(2p — 2) must be subtracted to account for 
blocks that satisfy multiple forms. Therefore 

/ ,\ (p — k)(p — k + 1) / \ /, , t2 P^~P\ / 1N 

ai+x,p{p -n + k) = ai+x,p{n) + [kp + k - k -\ — ) • ai+x,p{'^ + 1) 



k'^-k 
2 



ai+^,p(n + 2)-(2p-l)(2p-2). 



D 



This is equivalent to Theorem 5.10 of AUouche-Berthe [T], reproduced below in Theorem 2.5 





Blocks for k = 


Row 

mod 
P 





1 


2 




p-l 





61OOO. . .OO62OO 005„00 ...000 

O61OO. . .OOO62O 0006„0 ...000 

OO61O. . .000062 OOOOfen ...000 

0000 . . .61OOOO 6„_iOOOO. . .06nO 

0000 . . .O61OOO 06„_iOOO. . .006„ 


6n+l 










6n+l 















6n+l 




1 


61OOO . . .O6262OO 06„6„00 . . .006„+i 

6161OO. . .OO6262O 006„6„0 ... 000 

O6161O. . .0006262 0006„6n ... 000 

0000 . . . 61OOOO 6„_iOOOO ... 6„0 

0000 . . .6161OOO 6„_i6„_iOOO. . . 6„6„ 


6n+l 
6n+l 









6n+l 
6n+l 













6n+l 
6n+l 
















p-l 


6162(5262) {gibn+i){g?,bn+i){g2bn+i) 

(5261)6162 (55&n+l)(546„+l)(s'3^n+l) 

(53 61) (52 61) 61 {g(ibn+i){g^bn+i){gAbn+i) 

(526i)(53^i)(546i) (526„)6„6„+i 

61(5260(5361) (536„)(526n)6n 


6n+l 
(526n+l) 
(53fen+l) 

(526n+l) 
6n+l 


6n+2 

6n+l 

(526„+l) 

{g^K+i) 

(526„+i) 






(53&n+2) 
(546^+2) 
{gbK+2) 

bn+1 

(526„,+i) 



Table 1: Forms of blocks for the general case 1 + x with any prime p 



Theorem 2.5. For < fc < p — 1 and n > such that pn + k >3 



ai+x,2{pn + k + l) - ai+^_2(pn + k) ={p - k) [ai+^^2{n + 1) - ai+x,2{n) 

+ k{ai+x,2{n + 2) - ai+x,2{n + 1 



with starting points ai+a;,2(0) = 1, ai+a;,2(l) = P, ai+x,2(2) = p^ , and 01+2^,2(3) = - — ^ — - — . 



2.2.2 Accessible m-Blocks for c + x + x^ and prime p 



Table [21 provides ac+a;_|_^2 p(n) for small n and p. 



Using a method similar to the one we used for Theorem 2.4 the recursion relations appear to 
be those shown in Table |3j 
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Prime 


c 


a(n) 


2 


1 


2 4 8 4 25 36 53 70 92 114 


3 


1 


3 9 25 43 71 109 157 207 259 313 


3 


2 


3 9 25 61 105 165 233 321 417 533 


5 


1 


5 


25 121 393 673 929 1257 1761 2341 3097 


5 


2 


5 


25 125 393 689 953 1293 1801 2389 3145 


5 


3 


5 


25 117 385 657 905 1221 1713 2277 3017 


5 


4 


5 


25 101 169 253 353 509 721 989 1313 


7 


1 


7 49 331 1285 2137 2881 3859 



Table 2: a(n) for c + x + x^ 



p 


c 


Recursion 


k 


initial 


2 


1 


2a(n)+2a(n+l) 
a(n)+2a(n+l)+a(n+2) 


8 


1,2,4,8,14,25 


3 


1 


6a(n)+3a(n+l) 

3a(n)+6a(n+l) 

a(n)+7a(n+l)+a(n+2) 


20 


1,3,9,25 


3 


2 


4a(n)+4a(n+l)+a(n+2) 
2a(n)+5a(n+l)+2a(n+2) 
a(n)+4a(n+l)+4a(n+2) 


32 


1,3,9,25,61,105 


5 


1 


9a(n)+12a(n+l)+4a(n+2) 

6a(n)+13a(n+l)+6a(n+2) 

4a(n)+12a(n+l)+9a(n+2) 

2a(n) + 10a(n+l)+12a(n+2)+a(n+3) 

a(n)+12a(n+l)+10a(n+2)+2a(n+3) 


152 


1,5,25,121,393,673 


5 


2 


9a(n)+12a(n+l)+4a(n+2) 

6a(n)+13a(n+l)+6a(n+2) 

4a(n)+12a(n+l)+9a(n+2) 

2a(n) + 10a(n+l)+12a(n+2)+a(n+3) 

a(n)+12a(n+l)+10a(n+2)+2a(n+3) 


152 


1,5,25,125,393,689 


5 


3 


9a(n)+12a(n+l)+4a(n+2) 

6a(n)+13a(n+l)+6a(n+2) 

4a(n)+12a(n+l)+9a(n+2) 

2a(n) + 10a(n+l)+12a(n+2)+a(n+3) 

a(n)+12a(n+l)+10a(n+2)+2a(n+3) 


152 


1,5,25,117,385,657 


5 


4 


15a(n) + 10a(n+l) 

10a(n) + 15a(n+l) 

6a(n) + 18a(n+l)+a(n+2) 

3a(n)+19a(n+l)+3a(n+2) 

a(n)+18a(n+l)+6a(n+2) 


72 


1,5,25,101,169 



Table 3: Recursions for c + x + x^ 



We see that for p > 2, ac_,_2,_,_^2 p(n) = ai^x,p{'n) if c = | (mod p) because c+x + x^ = (l + x/2)^. 
Furthermore, we arrive at Conjecture |2.6[ 



Conjecture 2.6. For c ^ ^ (mod 5), the recursion for ai_^_r^_^_^2 Jn) is independent of c. Only the 
initial terms of the recursion depend on c. 

2.3 Closed form for a{n) 
Theorem 2.7. ai^,j:,2{m) = rn? — m + 2. 



Proof. Theorem 2.4 provides the recursion relation of ai^x,2{'2n) = 2>ai+x,2{n) + ai+x,2{n + 1) — 6 
and aiJ^x,2{n) = aiJ^x,2{n) + 'iaiJ^x,2{n + 1). We can find the starting points of 01+^.^2(1) = 2 and 
ai+a;,2(2) = 4 through inspection. This uniquely defines the sequence of accessible TTi-blocks. It is 
easy to show that the equation ai+a;(m) = m? — m + 2 satisfies both recursion relations through 
substitution, and also satisfies 01+^^2(1) = 2 and ai+x,2(2) =4. D 

This matches Remark 5.14 of [1]. 

2.3.1 Generating Functions for a(n) 

Using recursion relations, we can find the generating functions 9f{x),p ^^^ P > 3. 
Theorem 2.8. 



gi+x,p{z) = Y^ ai+x,p{^ 

l + (p-3)z + (/-3p + 3)z2 



n=0 

1 



(1-Z)3 



Proof. We have from Theorem 2.4 that for starting points a(0) = 1, a(l) = p, and a(2) = p the 



recursion relation is defined for pn + k > 2 as 



a{pn + k) ^ (P ^)(:P ''^^K {n) + {kp + k - k^ + ^^)a{n + 1) 
+ '^-^a{n + 2) - {2p - l){2p - 2). 



Adjusting for the k = 0,1 cases by replacing k with n + 2 gives 

a(pn + k + 2) ^ (p-fc-2Kp-fc-l) ^^^^ + (^^ _ 3^ _ ^2 _ 3 + ?!±i^)a(n + 1) 

(k + l)(k + 2) , , , 

+ ^ ^7 ^«(" + 2) - (2p - l)(2p - 2). 



To adjust for the case when p,k = 0, we define the recursion relation to have an additional term 
of '^~ 2 — ^^(0) + 2 — ^"^(1) ~ (2p — l)(2j' — 2) subtracted from the right hand side for only 

the case oi p,k = 0. 

We multiply through by 2;^"+^^ then sum over A; = to p— 1, then n = to cx). We also subtract 
from the right hand side of the sum the above mentioned additional term to account for the case 
oi p,k = 0. Defining h{x) = X) o(n + 2)z", we get 

n>0 

h{z) ={l + z + z' + ... + zP-'fhizP) + :-r^--. (p^zil - zf + 2/(1 _ ^)(4 _ 5^ + 2z^) 

2(1 - zy V 

+ 2(2 -2,z + 2,z^ -z^ - zP) - p{12 - 19z + 16^^ - 5z^ - Qz^ + 2z'^p] 
(2p-l)(2p-2) 

(I — zP)^ 1 

Therefore h(z) = -, 6rHz^) + Q{z) - (2p - l)(2p - 2) where 

(1 — z)-^ 1 — z 

Q{z) =-— ^— (p^z{l - zf + 2p2(i _ z){4 -5z + 2z^) + 2(2 - 3z + 3z^ - z^ - zP) 
2(1 - z)-^ \ 



- p{l2 - 19z + 16z^ - 5z^ - 6zP + 2z'^p) 
We then define u(z) = {l-zfh{z) and R{z) = (5(z)(l-z)^- (2p-l)(2p-2)(l-z)^. Iteratively 



substituting gives u{z) = u{zP°°) + J2 R{zp') = a{2) + J2 -R(^^'), or h{z) = TT^(a{2)+J2 R{z^' 

i>0 i>0 y ) \ j>o 

Note that 



i>0 i>0 



+ {p-2){p- ifz^ + 2(3p - l)zP - 2pz^P 

(3p - l)z - pz'^) + i^ZllL ^ (^pz'P'^ _ 2{p - l)z'^P* + {p- 2)z3p') . 



i>0 



Therefore 



g{z) =a(0) + a{l)z + 2;^/i(z) 

II I 2 i>0 

=1 +pz + z 



3 

2 

j>0 



(l-z)= 
1 + (p - 3)z + (p2 _ 3p + 3)^2 + ^2 (P:_i)! j2 (pzP' - 2{p - 1)z^p' + {p - 2)z^p') 



1-z) 



D 



Example 2.9. Setting p = 3 in Theorem 2.8 and noting that the z^P further reduces when p = 3 
provides 



-| (JXJ 

gi+x,3{z) = .^_ .3 (l + 3^2 - 2^3 + 8z2 ^(z3' - Z 



i=0 



Example 2.10. Setting p = 5 in Theorem 2.8 provides 



1 oo 

^- - ^".2 I Q^2Y^/c-^5' Q,2-5' I Q 3-5* 



5i+x.,5(2) = .^_ .3 (1 + 2z + 13^2 + 8^2 ^(5^5' - 8z2-5" + 3/- 



i=0 



We can use a similar proof to find further generating functions gx),p{z) from the recursion 
relations for af(^x),piji)- 

Theorem 2.11. 

l+2z3 + 2z5_^6+g(^2»_^3.2») 

10 



Based on the recursions in Table |3] and the method provided in Theorem 2.8, we arrive at 



Conjecture 2.12 which is confirmed for p = 3, 5. 



Conjecture 2.12. Fore ^ -r (mod p), the functional equation for the generating function g^. 



+x+x^ ,p 



{z) 



IS 



9c+x+x'^ ,p\Z) — 7~y9c-\-x+x'^,p\Z ) — Q{z) — — — , 

where r{z) = (1 — z^){l — z)"^ and Q{z) is some polynomial. 

Conjecture 2.13. For any f{x) and p, the generating function gfM.pi^) satisfies the equation 
r{z)gff^\p(z) = r(z'^)gff^\p(z'f) + b{z) for some polynomials r(z) and b{z) depending on f{x) and 
p. 



2.4 Limits of 



a[n] 



n^ 



Using the generating functions, we can find the asymptotic behavior of a{n) as n goes to infinity. 
Inspired by the quadratic nature of Theorem 



2.7 



we examine the behavior of 



a(n) 



Theorem 2.14. For f{x) = 1 + x and any prime p>2>, 



hm 

n— >oo 



ai+x,p{n) 



p^{p — 5){p — 1) 

2(p + l) V"' ' p{p-5) 



( , P+M ^ , (p-l)(p^-7p + 4) 
[x H — -, ^ ) + 



n^ 



+ 



-{p - l)(7p3 - 8p2 - 9p + 18) 

4(p + 1) 
(p-l)(p5 + 5p'^-8j 



2(:P-5) 

(p + l)(3p2-7p + 6)N^2 
7p3 _ 8p2 _ 9p + 18 



^ <x<\ 

p — — 3 



15p2 + 39p _ 18) 



2(7p3 - 8p2 - 9p + 18) 



(p-2)(p-l)(p^ + 2y + 5) / {p+l? 

4(p +1) V p2 _^ 2p + 5 

{p- l){p^ + Ap^ + ?,p - 4) 



+ 



2(p2 + 2p + 5) 



3 < X < 2 



\<X<1 



where n = \ — I and the limit as n —>■ oo is with constant x and k ^ oo. 

X 



Remark 2.15. The first polynomial from Theorem 2.14 corresponding to - < x < ^ should he 



understood in the sense of the limit for p = 5 as we divide by (p — 5). In this case the polynomial 
is not quadratic but actually the linear polynomial 20x + 8. 
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Proof. Theorem 2.8 states that 



9iz) = Yl ai+x,p{n)z 

n>0 
1 






(1-Z)3 



i>0 



Let y 6(n)z" = , ^ ,, y fpzP" - 2(p - 1)z^p' + (p-2}z 
— ' (1 — z) — ' ^ 

n>0 *^ '' i>0 

k 

Therefore, with the Hmit of n = [^J — ^ oo taken with fixed x and k — t- oo, we have 

n>0 ^^ ^> n>0 ^ 

Jrm a(n) = Jim (^^^^^ + ^^^^ + p + ^^6(n) j 

hm^=(^^^+(^^^lim^. 
Therefore, because they act similarly, we can find the asymptotics of "^^ by understanding the 



n^ 



behavior -^. We can rewrite Yl, b{n)z'^ as 
" n>0 



- / P^" (n-p-)(n-p--l) _ 2^^ _ ^^ '^" (n-2p')(n-2p'-l) 

n=0 ^ i=0 i=0 



+(j,-2) z ^"-'"■^^"-'"■-'K" 



2=0 



From this we see that 



i=0 \ ^ / i=0 V ^ 

+(?-2)g"( '"-^'''"7^'''-" ). 
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Therefore 



P'^" / ^i „J I 1 \ 2p'<n 



'fi|(a-^>-'^.)-(.-)i:(a-?.a-^.) 

fc 

Let n = [— J . We can neglect the 1 in the second factor (it creates a change that goes to zero as 
k —7- cx)), so we get 

j=0 1=0 j=0 

Note that if x [o, 1] then there is ?n, G Z such that p^x G [-, 1], so we can assume - < x < 1. 

fc 

Ignoring the floor for simplicity, we set n = — . Therefore we get 



X I i=0 j=0 i=0 



(?) 



When examining the upper limits of the three sums, we find that we therefore have 3 cases: 
i<x<i,i<x<^,^<x<l. For the first sum, p* < ^ gives i < /c + 1 for x = -^ and i < k 
for X = o, 2) 1- For the second sum, p* < f- gives i < A; for x = ^, o) s ^^'-^ i < A; — 1 for x = 1. For 



the third sum, p* < |^ gives i < k for x = -, o and i < k — 1 for x = 2, 1- Note that the limit is 

fc 

taken along the subsequences of the form [—J with fixed x and k — > oo. Also note that the limiting 
function does not change if x is replaced by p • x. 
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For the first case of - < x < g, we find that 



{ — ) ^ i=0 i=0 j=0 

2 1 1 

=(P - 5) „ ^^ x^ + 2 ^x 

p^ — 1 p — 1 

T K"-) T // ^\P ~^ 2 ^P ~ r?' 

iim — ^ = iim (p — 5) — ^ — ^^ — X + 2 ^— x 

n^oo n^ k^oo \ p^ — 1 p — 1 

ain) p'^ip- 5)(p-l) / p+1 \2 (p- l)(p'^ -7p + 4) 
lim — F- = — — xH — + 



ra-s-oo n^ 2(p + l) V p(p — 5)/ 2(p — 5) 

For the case of o < x < ^ we similarly find that because 



S^ = f E(i - p'-'-f - (^ - 1) E(i - v^^x)^ + "^ E(i - ^p'^'-f^ 

V^j j=0 1=0 i=0 



the limit of 

,. a(n) -(p-l)(7p^ -8p'^ -9p+18) f (p + l)(3p'^ - 7p + 6)\'^ (p-4)(p-lf 
rn-oo n? 4(p + 1) V 7p3 _ 8p2 _ 9p + 18 y 4 

(p+l)(p-l)(3p^-7p + 6)^ 
"^ 4(7p3 - 8p2 - 9p + 18) 

Similarly for the case of 2 < x < 1 we find that because 

fc-l o fe-i 



S^ = I E(i - p^-'^)' - (P - 1) E(i - v-^x)2 + "^ E(i - ^p'-'^)' 

(^) j=0 i=0 j=0 



one has 



a(n) (p-2)(p-l)(p'^ + 2p + b) / (p+1)^ \2 (p - l)(p3 + 4p2 + 3p - 4) 

n-^00 IT? 4(p + l) V |,2_^2p + 5/ 2(p2 + 2p + 5) 



D 
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Corollary 2.16. For the polynomial 1 + x and p>2>, 



lim inf °^+"'P(^^ = (p-l)(p' + V+3p-4) 

n-s>oo 7T,2 2{p^ + 2p + 5) 



lim sup 



ai+x,p(n) _{p- l){p^ + 5/ - 8p3 - I5p2 + 39^ _ is) 



n^ 



2(7p3-8p2-9p+l8) 



2.14 



is when x 



Proof. The maximum of Theorem 

p^ + 2p + l 

X = -^ . 

p^ + 2p + b 

We can also apply this to other aj(^)p(n). 
Theorem 2.17. For polynomial 1 + x + x"^ and prime 2, 



3p3 _ 4^2 _ p _|_ g 



and the minimum is when 



D 



-.. ai+^+a.2 2(n) 
lim ^ 



- + -X x^ 4 < X < 4 

4 2 12 2 - - 3 



n^ 



3 1 7 2 

12-4^+48^ 



< x< 1 



Furthermore, the upper and lower limits of ^+"^+2 c^^e - and — respectively. 



The proof of Theorem 2.17 is similar to the proof of Theorem 2.14 



Using the recursion relations, we computed the upper and lower limits of ^^"'•'2' — for sufficiently 
large n for several f{x) and p, The oscillatory nature of this sequence for large n stabilizing to a 
periodic function in log{x) is illustrated by Figure pi 




5.5 6.0 6.5 7.0 7.5 8.0 
(a) 1 -\- X (mod 3) 



4.0 4.5 5.0 5.5 

(b) 1 + a; (mod 5) 



(c) 1 + x + x^ (mod 2) 



Figure 3: — —^ with the x axis showing log„m, 
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This matches a prior result expressed in Lemma 5.15 by [T], which states that for large n, there 
exists constants ci and C2 such that cin? < a{n) < C2'n?- The limits given by Corollary 2.16 provide 
sharp values of ci and C2. [j 

3 Counting Coefficients 

3.1 Definitions 

For a polynomial /(x), prime p, and positive integer a < p — 1, we define qf(x),p{^^'^) to be the 
number of occurrences of a among the coefficients of f{x)^ reduced modulo p. Similarly, we define 
Qf(x),p{^) to be the total number of nonzero coefficients of f{x) . We then define rjiri.\p{n,a) = 

n—l n—1 

J2 Qf{x),p{h'^) and rj(^j.)^p(n) = J2 Qf(x),p{i)- We search for a quick method for calculating both 

i=0 i=0 

qf(x),p{k, a) and the asymptotic behavior of rf(^^-^,p{n, a) for large n. 

3.2 Willson Method 

Willson [B] describes an algorithm for computing the value of ^/{a;), 2(^)1 which is provided in 
Theorem 13. 1[ 

Theorem 3.1 (Willson's Method). For some polynomial f{x) with m,axim,um, degree d, there exists 
a matrix B, row vector u, and column vector v each of size 2°^ — 1 such that u- B^ -v = rf(x),2{'^^)- 

Amdeberhan-Stanley 0| describes a similar and related algorithm for calculating the number of 
each coefficient a for any power k for general f{x) and p, namely qf(x),p{k, a). Willson also analyzed 

the case of p > 2 in [7] . 

p 2" 
Example 3.2. For l + x + x"^ mod2, B = 1 1 2 . Note that the largest eigenvalue of this matrix 

[l 1 oj 

is 1 + \/5. 

Theorem 3.3. The matrix B is the sum of four matrices, each of which corresponds to a self- 
mapping of the set X = F2[x\/x'^ \ 0. 



^Strictly speaking, fore these sharp values, we may not have cin^ < a{n) < C2n^, but for any S > we have 
(ci — (5)n^ < a{n) < (c2 + S)n'^ for large enough n. 
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Theorem 3.3 follows easily from Willson [S]. 

Remark 3.4. The size of the matrix B can he made sm^aller only by using accessible blocks, as 
explained in Wilson [U^. 

3.3 Eigenvalue Analysis 

The matrix B has nonnegative entries and is irreducible. Following Willson pi, define A to be 
the Perron-Frobenius eigenvalue of B, i.e., the largest positive eigenvalue of B (it exists by the 
Perron- Frobenius theorem). We define A(/) to be the value of A for the polynomial /(x). We can 
approximate the value of 'ffM„{p'^,a) with A^ because the entries of B^ grow as a constant times 
A^ 

Example 3.5. For f{x) = 1 + x and p = 2, A = 3 because B = [3]. In this case A corresponds 
exactly to the scaling of the number of nonzero coefficients when doubling the number of rows, 
namely ri+a;,2(2A;) = 3 • ri+^^2{k). 

When examining the eigenvalues, we note that there are multiple transformations of a polyno- 
mial that does not change A. 

Theorem 3.6. We define the polynomials f{x) and g{x) to be similar if we can transform f{x) 
into g{x) through a combination of the transformations f{cx) and cf{x) with integer 1 < c < p, 
x^f{x) with integer c > 0, f{x'^) with integer c > 1, x -^^f' f{x^'^), and f{xY with integer c > 1. 
Any two similar polynomials have the same A. 

Proof. Because the transformations /(c • x), f{x^),x'^ ■ f{x), c ■ f{x), and flipping a polynomial do 
not change the number of nonzero coefficients of a polynomial, A do not change. Furthermore, 
because f{xY is every c*'* row, the ratios over the long term of the sums of total number of nonzero 
coefficients does not change, so A is the same. Namely, let Q/(x)(^) be the number of nonzero 
coefficients of /(x)". Therefore g/(x)(^ + 1) ^ C ■ <lf(x){n), where C is the number of nonzero 
coefficients of f{x). This means that 

k-n—l n—\ 

rfi.){k -71)= Y. qf(^,)U) < E(l + ^ + • • • + C^'"')9/W(i ■ k) < {1 + C + . . . + C7^-i)r^(,).(n). 
j=o i=o 
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Polynomial 


A 


d 


Polynomial 


A 


d 


l + x 


3 


1 


1 + X + x^ 


3.45686 


20 


l + x + x^ 


3.23607 


2 


1 + X + x^ + x^ 


3.49009 


20 


l + x + x^ 


3.31142 


4 


1 + X + x^ + x*^ 


3.50478 


10 


l + x + x"^ 


3.33159 


5 


1 + x^ + x^ + x'5 


3.53521 


20 


l + x + x"^ + x'^ 


3.3788 


7 


1 + X + X^ + X^ + x*^ 


3.53141 


19 


l + x + x'^ + x^ 


3.47662 


4 


1 + X + x'' + x^ 


3.50468 


17 


l + x + x"^ + x'^ + x^ 


3.45729 


4 


1 + X + x^ + x*^ + x*" 


3.55002 


19 


l + x + x^ 


3.35174 


10 


1 + X + x^ + x'^ + x*" 


3.59415 


16 


l + x"^ + x^ 


3.46127 


12 


1 + x^ + x^ + x'* + x"^ 


3.53665 


15 


l + x + x"^ + x^ 


3.49563 


7 


1 + X + x^ + x^ + x"^ + x*^ 


3.59043 


11 


l + x + x'^ + x^ 


3.45469 


12 


1 + X + x^ + x^ 


3.54536 


14 


1 + x^ + x^ + x^ 


3.46639 


5 


1 + X + x^ + x^ + x*" 


3.50809 


18 


1 + X + X^ + X^ + x^ 


3.5229 


14 


1 + X + x^ + x^ + x^ + x*^ 


3.57066 


17 


1 + X + x^ + x'' + x^ 


3.47168 


11 


1 + X + x^ + x"^ + x^ + x*^ 


3.49995 


6 


1 + X + X^ + X^ + X'* + x^ 


3.52951 


6 


1 + X + X^ + X^ + X^ + x^ + x^ 


3.5598 


6 



Table 4: A and the degree of its minimal polynomial for p = 2 and deg(/(x)) < 6 



This implies that A(/) < \{f^). Similarly since qf(x)U • k — i) > C ^qf(^x)U ' k), we can show that 



\if) < A(/). Therefore A(/) = X{f). 



a 



3.3.1 Values of A where p = 2 

We calculate A for polynomials with p = 2. We also find the minimal polynomial of A. Provided 
are A and the degree d of its minimal polynomial for non-similar polynomials with degree of up to 
6, although we had calculated for deg{f) < 9. 

We see that A is between 3 and 4. We form several conjectures on the bounds of A. 

Conjecture 3.7. When p = 2, A > 3. Furthermore, A = 3 only for polynomials similar to l + x. // 
p = 2 and A > 3, then A > 1 + \/5- Furthermore, A = 1 + \/5 only if /(x) is similar to 1 + x + x^. 

Question 3.8. Is it true that A(/) = X{g) if and only if f{x) and g{x) are similar in terms of the 
transformations described in Theorem\3.6\? 



Theorem 3.9. For some polynomial /(x) with degree at most 2^ and p = 2, 

A(/)<4(l-2^)^. 
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Proof. Define k such that the degree of f{x) is at most 2'^, with p = 2. From Theorem 3.3 we can 



draw an oriented graph whose vertices are elements of X and whose edges correspond to the four 

maps. Therefore there are exactly four edges coming out of each vertex. Therefore if Q{n) is the 

los (D(ti) 
number of paths in the graph of length n, we have log A = limsup . From the definition 

n— >oo n 



of Willson's method, Theorem 3.1 two of the four mappings correspond to g{x) -^ g{x'^) and 
g{x) — ^ X • (?(x^). Assume deg(/(x)) = 2 . Then a path starting from any g{x) and moving first to 
X ■ g{x'^) then alternating in any way between the two mappings leads to after k + 1 steps. So the 
number of such paths of length /c + 1 is 2 . So the number of paths of length k + 1 from any point 
that avoids is at most 4^"^ — 2 . Thus the number of such paths of length n • (A; + 1) is at most 
(4^+1 _ 2^=)". This gives us the bound of A < 4(1 - ^)'^ . D 

For k = 0, the only polynomial is 1 + x, so the bound A < 4(1 — |)^ = 3 is sharp. However, for 
k = 1 the bound tells us that A < \/l4 which is not sharp. Furthermore, this bound approaches 4 
as k approaches oo. 

Conjecture 3.10. Let A^ be the maximal A(/) for deg f < k. Then limi^^oo-^k = 4. 

Remark 3.11. Similarly for p > 2, one may conjecture that linik^oo-^k = P^ ■ 

Through computer analysis of A for p = 2 and deg{f{x)\ < 9, Conjecture 



3.12 



arises. 



Conjecture 3.12. The degree of the minimal polynomial of X is less than or equal to 2'^'^^^^' ^ for 
p = 2. 

4 Conclusion and Directions of Future Research 

Natural goals for further study of the phenomena examined in this paper include the following: 

• Obtain recursion relations, generating functions, and limiting functions as in Section 2 for 
aj(^)p(n) in the case deg(f{x)j > 1; 



Prove Conjecture 2.13 on the functional equation for the generating function for a/(a;),p(^); 
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Prove the conjectures in section 3 on the behavior of the eigenvalues A and obtain better 
upper bounds; 



Find, tighten, and explore the upper bound mentioned in Conjecture 3.12 



• Study the algebras generated by the four transformations composing the Willson matrices 
and find analogs for larger p. 
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